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SOLVING GAUGE-FIELD THEORY 
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The canonical front form Hamiltonian for non-Abelian SU(N) gauge theory in 3+1 
dimensions is mapped non-perturbatively on an effective Hamiltonian which acts 
only in the Fock space of a quark and an antiquark. The approach is based on 
the novel method of iterated resolvents and on discretized light-cone quantization, 
driven to the continuum limit. It is free of the usual Tamm-Dancoff truncations of 
the Fock space, rather the perturbative series are consistently resumed to all orders 
in the coupling constant. Emphasis is put on dealing with the many-body aspects of 
gauge field theory. The effective interaction turns out to be the kernel of an integral 
equation in the momentum space of a single quark, which is frame-independent 
and solvable on comparatively small computers. Important is that the higher 
Fock-space amplitudes can be retrieved self-consistently from these solutions. 


1 The Hamiltonian matrix 

In canonical field theory the four components of the energy-momentum vector 
P v commute and are constants of the motion. In the front form of Hamiltonian 
dynamics EJ they are denoted a by P v = (P + , P±,P~). Its spatial components 
P_l and P + are independent of the interaction and diagonal in momentum 
representation. Their eigenvalues are the sums of the single particle momenta, 
P + = Ylv + = and P± = Y^P-L- Each single particle has four-momentum 
= (p + ,p±,p~) and sits on its mass-shell p^P/j = m 2 . Each particle state 
“q” is then characterized by six quantum numbers q = (p + ,p± 1 A; c, /): the 
three spatial momenta, helicity, color and flavor. The temporal component 
P~ = 2Pk. depends on the interaction and is a very complicated non-diagonal 
operator0. It propagates the system in the light-cone time x + = x° + x 3 . The 
contraction of P M is the operator of invariant mass-squared, 

P M P M = P+P- - P 2 = H hC = H . (1) 

It is Lorentz invariant and referred toi somewhat improperly but conveniently 
as the ‘light-cone HamiLtonian’ Plc> or shortly H. One seeksaurepresentation 
in which H is diagonal Ena, H |\I>) = P|'I'). Introductory textstffi are available. 

By convenience, the (light-cone) Hamiltonian is split into four terms: H = 
T+V+F+S. The kinetic energy T survives the limit of the coupling constant g 
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Figure 1: The Hamiltonian matrix for a SU(N)-meson and a harmonic resolution K = 4. 
Only vertex diagrams are included. Zero matrices are marked by (■). 


going to zero. Since it is diagonal in Fock-space representation, its eigenvalue 
is the free invariant mass squared of the particular Fock state. The vertex 
interaction V is the relativistic interaction per se. It is linear in g and changes 
the particle number by 1. Those matrix elements of V which change the 
particle number by 3 are present in the instant form, but vanish for light-cone 
kinematics: The vacuum does not fluctuate. The instantaneous interactions 
F and S are consequences of working in the light-cone gauge A + = 0. They 
are proportional to g 2 . By definition, the seagull interaction S conserves the 
particle number and the fork interaction F changes it by 2. 

The eigenfunctions |\I/) are complicated superpositions of Fock-space pro¬ 
jections with contributions like (gg|\I/) or (qqg\'$>). The eigenvalue equation 
stands thus for an infinite set of coupled integral equations which are even dif¬ 
ficult to write down. However, if one works with periodic boundary conditions, 
with discretized light-cone quantization □ (DLCQ), the infinite set of coupled 
integral equations is mapped onto a finite set of coupled matrix equations, 
as illustrated in Fig. [I], As consequence of discretization, the Fock states are 
denumerable and orthonormal. The momentum fraction carried by the n-th 
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particle is x n = p^/P + . The sum of all momentum fractions and the sum of 
all transversal momenta is constrained to x n — 1 and Y^ n kn± = 0, with 
the intrinsic transversal momenta k± defined by p n± = x n P± + k n± . Since P + 
has only positive eigenvalues and since each particle has a lowest possible value 
of p + , the number of paRicles in a Fock state is limited for any fixed value of 
the harmonic resolutions .I K. 

Here then is the problem, the bottle neck of any Hamiltonian approach in 
field theory: The dimension of the Hamiltonian matrix increases exponentially 
fast, even if one regulates the transversal momenta by some suitable cut-off A, 
like the invariant-mass cut-off of Lepage and Brodsky □ or similar constructs. 
Suppose, the regularization procedure allows for 10 discrete momentum states 
in each direction. A single particle has then about 10 3 degrees of freedom. 
A Fock-space sector with n particles has then roughly 10 ra_1 different Fock 
states. Sector 13 with its 8 particles has thus about 10 21 and sector 1 (qq) 
about 10 3 Fock states. Deriving an effective interaction can be understood as 
reducing the dimension in a matrix diagonalization problem from 10 21 to 10 3 ! 

2 The effective interaction 

Effective interactions_are a well known tool in many-body physics. In field 
theory the method 1913 is known as the Tamm-Dancoff-approach. The Fock- 
space sectors |f) appear in a most natural way and the Hamiltonian matrix 
can be understood as a matrix of block matrices, whose rows and columns are 
enumerated by i = 1, 2,... IV like in Fig. [[]. The eigenvalue equation can then 
be written as a coupled block matrix equation: 

N 

= E (nl’F) for all* = 1,2,..., IV . (2) 

l=i 

The rows and columns can always be split into the P- space and the rest, the 
Q-space. The block matrix equation becomes 

(P\H\P) (P\V) + {P\H\Q) (Q\V) = E(P\V), (3) 

and (Q\H\P) (P\V) + (Q\H\Q) (Q\ty = E (Q|<F) . (4) 

Rewrite the second equation as (Q\E — H\Q) (Q|’F) = (Q\H\P) (P|^), and 
observe that the quadratic matrix (Q\E—H\Q) could be inverted to express the 
Q-space wave-function (QI'F) in terms of (P|^). The eigenvalue E is unknown 
at this point. One introduces therefore a redundant parameter ui, and defines 
Gq(uj) = [{Q\u} — -ff |Q}] -1 . The so obtained effective interaction 

(P\H eS (u)\P) = (P\H\P) + (P\H\Q) G q {lo) ( Q\H\P) . (5) 
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acts only in the P- space: H e ff(co) |$fc(w)) = Ek{to) |4>fc(u;)). Varying lo one 
generates a set of energy functions Ek(to). Whenever one finds a solution to 
the fix-point equation Ek(co) = w, one has found one of the true eigenvalues 
and eigenfunctions. The advantage of working with an effective interaction is 
of analytical nature to the extent that resolvents can be approximated system¬ 
atically. The two resolvents 

“ d G ° M = (Q|u- t \ q ) ’ (6) 

defined once with and once without the non-diagonal interaction in the Hamil¬ 
tonian H = T + U, respectively, are identically related by Gq{lo) = Gq(co) + 
Gq(u) U Gq(lo ), or by the infinite series of perturbation theory Gq{lo) = 
G 0 (lo) + G 0 (ui)UGq(lo) + G 0 (lo)UG 0 (u>)UG 0 (lo) + ... . The point is, of course, 
that the kinetic energy T is a diagonal operator which can be trivially inverted 
to get the unperturbed resolvent Gq(u>). 

Albeit exact in principle, the Tamm-Dancoff-approach (TDA) suffers in 
practice from two serious defects: (1) The approach is technically useful only, 
if one truncates the perturbative series to its very first term. This destroys 
Lorentz and gauge invariance. (2) If one identifies lo with the eigenvalue, as 
one should, the effective [injteraction has a non-integrable singularity, both in 
the original instant formBQ and in the front forrnBB. 

The above can be understood as to reduce the block matrix dimension 
from 2 to 1. But there is no need to identify the P -space with the lowest 
sector. One also can choose the Q-space identical with last sector and the 
P -space with the rest, P = 1 — Q: The same steps as above reduce then the 
block matrix dimension from N to N — 1. The effective interaction acts in 
the now smaller space. This procedure can of course be iterated, but tlpn 
one has to deal with ‘resolvents of resolvents’, or with iterated resolvents liB. 
Ultimately, one arrives at block matrix dimension 1 where the procedure stops: 
The effective interaction in the Fock-space sector with only one quark and one 
antiquark is defined unambiguously. 

Suppose, in the course of this reduction, one has arrived at block matrix 
dimension n, with 1 < n < N. Denote the corresponding effective interaction 
H n (co). The eigenvalue problem reads then 



^2(i\H n (uj)\j)(j\'fr(u})) = E{w) (il’F(w)} , for i = 1,2,..., n. (7) 

j =i 

Observe that i and j refer here to sector numbers. Now, like in the above, 
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define the resolvent of the sector Hamiltonian H n (u>) by 


Gn{w) — -j —r- v f \I \ 

( n\oj — H n {u>)\n) 

n —1 

thus (nl^o;)) = G n {ui) J2( n \ H n(u;)\j) 0 'l*M> 

i=1 

The effective interaction in the (n — l)-space becomes then 

H n -iH = H n (u) + H n (u)G n (w)H n (w) (10) 

for every block matrix element (i\H n -i(uj)\j). Everything proceeds like in 
above, including the fixed point equation E(u) = to. But one has achieved 
much more: Eq.(^|) is a recursion relation which holds for all 1 < n < N\ 
Since one has started from the bare Hamiltonian in the last sector, one has to 
convene that Hn = H. The rest is algebra and interpretation. 

Applying the method to the block matrix structure of QCD, as displayed 
in Fig. [I], is particularly easy and transparent. By definition, the last sector 
contains only the diagonal kinetic energy, thus H 13 = T 13 . Its resolvent is 
calculated trivially. Then H 12 can be constructed unambiguously, followed by 
-Hn, and so on, until one arrives at sector 1. Grouping the so obtained results 
in a different order, one finds for the sectors with one gg-pair: 


H qg = 


= 2 i + VG 3 V + VG 3 VG 2 VG 3 V , 

( 11 ) 

Hqqg = 

h 3 

= T 3 + VG 6 V+ VG 6 VG 5 VG e V+ VG 4 V , 

( 12 ) 

Hqqgg = 

He 

= T 6 + VG10V + VG10VG9VG10V + VG 7 V. 

(13) 


The quark-gluon content of the respective sectors is added here for an easier 
identification. Correspondingly, one obtains for the sectors with two qq -pairs 

H qm = H 4 = T 4 + VG 7 V + VG 7 VG 6 VG 7 V , (14) 

H qmg = H 7 = Tr + VGnV + VGuVGwVGnV + VGsV.ilb) 

In the pure glue sectors, the structure is even simpler: 

Hgg = H 2 = T 2 + VG 3 V + VG 5 V , (16) 

H ggg = H 5 = T 5 + VG 6 V + VG g V . (17) 

Note that the above relations are exact. They liold 0 for an arbitrarily large K 
and thus in the continuum limit. Note also that the vertex interaction appears 
only in even pairs, typically in the combination VGV. This is the deeper 


(8) 

(9) 


5 



■ 




3 

2 

u i,o 

■ 

3 



1 3 

U 1,1 

3 1 

U 1.2 


Figure 2: The three graphs of the ef¬ 

fective interaction in the qq- space. — The 
lower two graphs correspond to the chain 
U = VG 3 V, the upper corresponds to U a = 
VG 3 VG 2 VG 3 V. Propagator boxes are repre¬ 
sented by vertical dashed lines, with the sub¬ 
script ‘n’ referring to the sector numbers. 


effective Hamiltonian H e 3 = H qq as 
As a net result the interaction scatters a quark with helicity \ q 


reason why the instantaneous interactions can be implemented^ ex post. The 

given by Eq. (0) is illustrated in Fig. § 

and four- 

momentum p = (xP + , xP± + k±,p~) into a state with A^ and four-momentum 
p' = {x'P+, x , P ± +k^p'-). 

In the continuum limit, the resolvents are replaced by propagators and the 
eigenvalue problem H e s\ip) = M 2 \ip) becomes an integral equation 


Mb (X, /u]_? A q 1 Agl^fr) — 


m a 


■kl 


m g 


■kl 


1 — x 


(x, k±; Xq, Xq\l/jb) 


52 [ dx'd 2 kl(x,k±;X q ,Xq\U+ U a \x',kl-,X' q ,X'q) {x',kl-,X' q ,X'q\ipb) ( 18 ) 


The domain D restricts integration in line with regularization. The effective 
potential U is diagrammatically defined in diagram of Fig. || and cannot 
change the flavor of the quark. It includes all fine and hyperfine interactions. 
Diagram U\ ^ contributes to the effective quark mass to. The flavor-changing 
annihilation interaction U a = U\ t o is probably of less importance in a first 
assault. The eigenvalues M 2 are the invariant mass 2 of a physical particle 
and the corresponding wave-function (x, k±; X q , Xq\4>b) gives the probability 
amplitudes for finding a flavored quark with momentum fraction x, intrinsic 
transverse momentum k± and helicity X q . These are boost-invariant quantities. 
The eigenfunctions represent the normalized projections of the full eigenfunc¬ 
tion |\E') onto the Fock states | q\q) = b^ q d\\vac). 

The knowledge of ipb is sufficient to retrieve all desired Fock-space com¬ 
ponents of the total wave-function. The key is the upwards recursion relation 
Eq.(||). Obviously, one can express the higher Fock-space components (ra|\I/) 
as functionals of il> qq by a finite series of quadratures, i.e. of matrix multi¬ 
plications or of momentum-space integrations. One need not solve another 
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eigenvalue problem. I show this by calculating the probability amplitude for 
a | gg)- or a \qqg)-st&te in an eigenstate of the full Hamiltonian. The first two 
equations of the recursive set in Eq.(|9j) are 

<2|*> = G 2 (2|tf 2 |l)(l|T), (19) 

and (3|tf) = G 3 (3|tf 3 |l)(l|^)+G 3 (3|fi 3 |2}<2|'F)- (20) 

The sector Hamiltonians H n have to be substituted from Eqs. and (p~6|). 
In taking block matrix elements of them, the formal expressions are simplified 
considerably since many of the Hamiltonian blocks in Fig. |l] are zero. One 
thus gets simply (2|j? 2 |l) = (2|V'G 3 V’|1) and therefore (2|*F) = G 2 VG 3 V { l|'k}. 
Substituting this into Eq.(jl^) gives (3|T) = G 3 F(1|'F) + G 3 V'G 2 VG 3 V'(l|\I r ). 
These findings can be summarized more succinctly as 


and 


I ^gg) 
I ^qqg) 


G gg VG qqg V \lfqq)i (21) 

GqqgV \lfrqq) + G qq g V G gg V G qq g V \lpqq) ■ (22) 


Note that the above relations are exact. The finite number of terms is in 
strong contrast to the infinite number of terms in perturbative series. Iterated 
resolvents sum the perturbative series to all orders in closed form. 

3 The mean field and the vertex functions 

Let us discuss in some greater detail the structure of the sector Hamiltoni¬ 
ans, particularly of H qqg in Eq.(n2|). The corresponding graphs are displayed 
diagrammatically in Figs. [| and]J Those in Fig. |] differ from look those in 
Fig. |2j only by an additional gluon. The gluon does not change quantum num¬ 
bers under impact of the interaction and acts like a spectator. The graphs in 
Fig. |] will be referred to as the ‘spectator interaction’ U 3 . Correspondingly, 
the graphs of Fig. ^will be referred to as the ‘participant interaction’ C/ 3 . The 
gluon quantum numbers are changed and the gluon participates in the inter¬ 
action. The analogous separation into spectator and participant interactions 
can be made in all quark-pair-glue sector Hamiltonians: 

H n =T n + Un + Un, for n = 3,6,10,15,... . (23) 

Note that the spectator interaction has the same diagrams as in Fig. except 
for the additional free and non-interacting gluons. 

The essence is that the spectator interactions U n define mean fields, and 
that they can be associated with their own resolvents G n : 


G n = 


u - T n - U n 


(24) 
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Figure 3: The three graphs of the spectator in- 

teraction in the 95 9 -space. Note the role of the Fi S ure 4: Some six S ra P hs of the P ar ' 
gluon as a spectator. ticipant interaction in the qqg- space. 


while G n was 

G n = -. (25) 

u-T n -U n -U n 

Both are related by G n = G n + Gn U n G n . Equivalently, they can be written 
as an infinite series G n = G n -\-G n U n G n -\-G n U n G n U n G n -\- ... , as usual. The 
main difference to the above conventional series is, that there the ‘unperturbed 
propagator’ Gq(lu) refers to the system without interactions while here the 
‘unperturbed propagators’ G n include the mean fields U n . One therefore deals 
here with ‘perturbation theory in medium’. Note that the present series differ 
from the above also with respect to the physics: The system stays in sector n. 
This allows for the identical rearrangement 


G 


n 


with R n 


Rn G n R n , 

1 + ^G n U n + —G n U n G n U n + ... 

z o 

1 


1 — G n U n 


(26) 


(27) 


This can be verified order by order and has, to our recollection, not been given 
before. The square matrices R n are always sandwiched between a quark-pair- 
glue resolvent G and the vertex V. It is therefore meaningful to introduce the 
effective vertex V = R n V or V = R n V, thus 

VG n V = VG n V . 
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(28) 





































The vertex operators R n come with same power as the coupling constant. 
They are a genuine part of the interaction and accumulate the many-body 
aspects manifesting themselves as non-point-like vertices. Finally, rewriting 
systematically Eqs. ©-©> one gets 


H 6 

= T 6 + VG w V + VG w VG 9 VG 10 V , 

(29) 

h 3 

= t 3 + vg 6 v+ vg 6 vg 5 vg 6 v , 

(30) 

H 1 = H 1 

= Ti + V G 3 V + V G3VG2V G 3 V . 

(31) 


Instead of being similar, the quark-pair-glue sector Hamiltonians H n = T n +U n 
now apparently all look the same. Note that the resolvents G2 , G 5 and Gg carry 
no bar. They correspond to the pure glue sectors. The distinction between 
participants and spectators makes no sense there. 


4 The approximations 


Thus far the approach is formally exact. For calculating the effective interac¬ 
tion the two resolvents G 3 and G 2 and the vertex function R 3 are needed. How 
can they be approximated? Suppose to have solved the eigenvalue problem in 
the qq- space, 


'^2(q;q\H q g(uj)\q l - 1 q , )i>b{^)) = M%(u)(q;q\ip b (w)). (32) 


The eigenvalues are enumerated by b = 1,2,.... The corresponding eigenfunc¬ 
tions (q; q\il>b{u)) are a complete set. Despite working in the continuum limit, 
we continue to use summation symbols for the sake of a more compact notation. 
Suppose further that an 10 was found which has the same value as the lowest 
eigenvalue M 2 = Mf. The substitution u = M 2 will hence forward be done 
without explicitly mentioning. Next, ask for the eigenvalues and eigenfunctions 
in the qq g- space. One need not to solve another eigenvalue problem, since the 
result is known ahead of time! By construction, the gluon moves relative to 
the meson subject to momentum conservation. The eigenfunction is a thus a 
product state |V’b.s) = | ipb) < 8 > |<p s )- Parameterizing the gluons four-momentum 
as = ( yP + ,yP± + qj_,Pg ), the eigenvalues are 


M 2 _ + g , g 

“ (1 - V) + V ■ 


(33) 


Knowing the eigenvalues and eigenfunctions, one can calculate the exact resol¬ 
vent. After a few identical rewritings one gets 


{q;q;g\G 3 \q ;q ;g) 


G3(g; q ; g) [(q; q; g\q ; <?'; g) - (q ; 9 ; g\A\q; q\ g 1 )] . (34) 
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The operator A can never become a Dirac-<5 function (q; q; g\q'; q'; g'), since 


A _ I / \ — M 2 ) ! , ! 

M * b V +y(M?~ ill 2 ) WU ' 


(35) 


It is therefore dropped, as an approximation. The proportionality coefficient 
is 

G 3 (q;q-,g) = ^ > with Q 2 = y 2 ^ 2 + q±- 

The same approximation yields for the two gluon propagator 

1 


(36) 


92 ) — 


M 2 — M 2 b ' 


(37) 


It is parameterized in terms of the glue ball mass M 2 b . 

For calculating the operator i?§, one addresses first to calculate the oper¬ 
ator B := —G^iyG 4 V + VGqV ). An inspection of Fig. |] yields that the loop 
diagrams (34) and (3fi) are the leading terms, by far. The have been calculated 
long ago by Thorn E3 for n c = 3 colors and massless quarks, using the per¬ 
turbative propagators. Repeating this calculations with the non-perturbative 
propagators G 4 and G§, which can be approximated in line with Eq.(j34|), one 
gets a diagonal B with 


B = --a s co In 


1 + 


1 A 2 ' 


(38) 


with a s = g 2 /47r and cq = (33 — 2iV/)/67t_.The latter number is not unfamil¬ 
iar from the work of Gross and WilczekEJ. Details will be given elsewhere. 
Applying the same procedure to even higher spaces gives again and again the 
same type of vacuum polarization diagrams, and thus a diagonal 

R 2 = -^-. (39) 

1 + - T5 - 


This continued fraction can be resumed to all orders by R 2 
and solved 


R 


2 

3 


2 


/ 


I+ 1 /I + a s Co In 

1 , 1 A 2 

|_ i+ (i-y) W\ 


Now, one can calculate analytically the effective interaction! 


1/(1 +BRI) 

(40) 
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5 Outlook 


The strength of the interaction (i? 3 ) tends to zero for increasing cut-off A —> oo. 
One may (or may not) interpret this result as asymptotic freedom. But this 
would be premature. The above is only the regulated effective interaction , al¬ 
beit being well defined for any fixed value of the cut-off A. The dependence 
on A must be removed by a future renormalization group analysis Eia. It looks 
as if the present explicit form-is almost ideally suited for an application of the 
Hamiltonian flow equations 1U. Finally, one should emphasize that the above 
approach is based on comparatively weak approximations being thoroughly 
non-perturbative: No smallness assumption on the coupling constant was ever 
needed, nor was the Fock space ever truncated. Gauge invariance is main¬ 
tained, and Lorentz covariance (including rotations) is strictly observed in the 
continuum limit. 
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